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Abstract 

In this work gelatinized potato starch is shown to retain the memory of past loading 
history. It exhibits a visco-elastic response which does not depend solely on instan¬ 
taneous conditions. A simple squeeze flow experiment is performed, where loading 
is done in two steps with a time lag r ~ seconds between the steps. The effect on 
the strain, of varying r is reproduced by a three element visco-elastic solid model. 
Non-linearity is introduced through a generalized calculus approach by incorporat¬ 
ing a non-integer order time derivative in the viscosity equation. A strain hardening 
proportional to the time lag between the two loading steps is also incorporated. This 
model reproduces the three salient features observed in the experiment, namely - 
the memory effect, slight initial oscillations in the strain as well as the long-time 
solid-like response. Dynamic visco-elasticity of the sample is also reported. 
Keywords: Spreading, Gelatinized starch. Visco-elasticity, Generalized calculus 


1 Introduction 


Complex fluids such as starch gels generally behave differently from both Newtonian 
fluids and Hookean elastic solids as discussed by Evans and Haisman(1980), Che et 
al. (2008) and Schoheld and Blair(1937). Their rheology involves visco-elastic or 
visco-plastic behavior. Generalized calculus is a useful technique for describing and 
analyzing such materials, as shown by Das(2011), Heymans and Bauwens( 1994)and 
Dutta Choudhury et al.(2012). One typical signature of their complexity is that in 
a stress-strain experiment their response at a certain time instant does not depend 
only on the stimulus acting at that instant, they ‘remember’ the history of previous 
treatment and respond accordingly. We describe a simple experiment which clearly 
demonstrates the effect of different histories of loading on a drop of potato starch 
squeezed between two parallel glass plates. 

Very simple squeeze-flow experiments on different fluids, where a drop of fluid 
in a Hele-Shaw cell is subjected to a load, have shown interesting results[7]. The 
experiment consists of placing a weight of the order of kilograms on a drop of fluid 
sandwiched between two transparent plates. Spreading of the fluid is measured by 
video-recording the area of contact between the plate and fluid. Fluids with non- 
Newtonian character show novel features like oscillations in the area of contact [6]. 
In the present experiment we report the effect of step-loading on a gelatinized starch 
solution. We show that when two loads are placed one after the other with a certain 
time lag, the asymptotic area of spreading is different from what is observed when 
the same two loads are applied simultaneously. We show that a three-element visco- 


elastic model, incorporating a non-Newtonian dastipot and a strain hardening of the 
elastic component, is successful in reproducing the experimental observations. The 
strain hardening is proportional to the time lag between the two successive loading 
steps. Dynamic visco-elasticity measurements have been done using a standard set 
up for supporting characterization of the sample. All experiments show that the 
sample exhibits long-time solid-like behaviour. 

2 Materials and Method 

2.1 Sample preparation 

The non-Newtonian fluid under study is an aqueous gel of potato starch (Loba- 
chemie, Mumbai). In our experiment, 2.5 g of potato starch is dissolved in 100 ml of 
distilled water and it is heated up for 10 minutes and boiled for 2 minutes, stirring 
continuously, allowing it to gelatinize. The solution is cooled for 1 hour and a pinch 
of dye is added to enhance the contrast. 

2.2 Squeeze-flow set up. 

A droplet of the non-Newtonian fluid is placed on a smooth glass plate of thickness 

1.3 cm and diameter 15.5 cm. The mass of the droplet varies between 0.04 gm to 
0.06 gm. A similar glass plate is placed on top of the drop. The weight of the upper 
plate is 0.56 kg. The entire experiment consists of three parts: 


Experiment 1: Pull loading 


30 s after placing the upper plate, a weight W is placed on it and a CCD 
camera (WATEC-202D, Japan) below the lower plate, records the spreading 
of the area of contact. The video recording is analysed using the software 
Image-Pro-Plus. 

• Experiment 2: Step-loading-A 

The procedure is repeated with a drop of equal volume, this time the loading 
is done in 2 steps. 30 s after placing the upper plate, weight Wi is placed on it. 
After a further interval of r = 5.3s another weight W 2 is added to the load on 
the upper plate. Here the hnal load W = Wi + W 2 . The strain development 
in this case for f > r is compared with the case for one-time loading. 

• Experiment 3: Step-loading-B 

The procedure for experiment 3 is repeated exactly, only the time lag r between 
placing Wi and W 2 on the upper plate is r = 12.3s. 

We now compare results of the above three experiments. Taking the instant of 
placing the hrst weight as f = 0, the fractional change in the area of contact of fluid 
and glass at time t 


e{t) = [A{t) - A(0)]M(0) (1) 

is measured using Image-Pro-Plus software. Here A{t) is the area at time t and 
A(0) the initial area. We identify this as the strain in the spreading experiment. It 
is difficult to make the initial area A(0) exactly equal in the experiments comparing 


different loading histories. To overcome this problem the experiment is repeated 
many times and only sequences where A(0) varies within 3% are chosen for analysis. 

2.3 Measurement of complex visco-elastic moduli 

Response of the sample to a sinusoidally varying dynamic load has been done for 
additional characterization of the visco-elastic behaviour. Measurements have been 
carried out at the Pharmaceutical Engineering Department, Jadavpur University 
using MCR102 SN81260812. 

3 Results 

3.1 Squeeze flow results 

The strain e{t) is plotted as a function of time for all 3 experiments. The results are 
summarized in hgureffTl). Two sets of data for e{t) versus t for each of experiments 
(1-3) are plotted to show the reproducibility of the results. The two sets for identical 
conditions are designated as data sets a and b. 

In Experiment 1, fU = 4.0 kg and in Experiments 2 and 3, Wi = W 2 = 2.0 kg. 
The time lag r between placing the two weights is 5.3 s in Experiment 2 and 12.3 s in 
Experiment 3. Figure([T]) clearly shows the deviation from both ideal Hookean elastic 
behavior as well as from the behaviour of an ideal Newtonian fluid. The strain does 
not reach the hnal value instantaneously as expected for a perfectly elastic material, 
neither does it rise from zero and continue to creep indehnitely like a perfect fluid. 
Moreover there is clear evidence of history dependence in strain for e{t) when t > r. 



Figure 1: Results for the strain variation with time in single step and two-step 
loading. Two data sets for each of the experiments are shown - solid circle and box 
are for single loading with 4kg, triangle up and diamond are for step loading with r 
= 5.3 s; triangle right and open circle are for step loading with r = 12.3 s. The solid 
line represents the calculated curves for 4 kg (single step loading), dashed line and 
broken line are for two-step loading with r = 5.3 s and 12.3 s respectively, using the 


3-parameter solid model. 









The final asymptotic strain is different in all three cases - being maximum when 
the loading is done at once and decreasing as r increases. Clearly, the continued 
loading changes the elastic and/or viscous properties of the material. We analyse 
the strain build-up using different visco-elastic models and see which model is best 
suited to reproduce experimental results. We try to understand the characteristics 
of the material through these results. 

3.2 Results for dynamical loading: the complex modulus 

The complex visco-elastic modulus is written as 

G\u) = G'{u) + iG"{u) (2) 

In the dynamic loading experiment the sample is subjected to a sinusoidally varying 
strain with frequency u rad/s and constant amplitude. Measuring the stress response 
which has a phase lag with the strain, gives the complex modulus G*. The real part 
G' gives the elastic contribution while the imaginary part G” represents the loss 
modulus, i.e. the viscous part. Graphs for G' and G” in the frequency range 0.1 
to 100 Hz are shown in figure (j2l). G' and G” both increase with frequency and the 
curve for G' always lies above G", showing a predominantly solid-like response. This 
type of behaviour is typical of polymer gelsjH |9] and starches [TO]. 



Figure 2: Experimentally measured storage and loss modulus. Circles with solid 
line (black) represent G' and squares with red solid line represent G". 







4 Visco-elastic modelling 


Modelling visco-elastic materials is usually done using energy conserving elements 
- like a spring and dissipative elements like a dash-pot. Different combinations of 
these basic elements lead to different constitutive stress-strain relations. 

The simplest visco-elastic models are the Maxwell model, where a spring is in 
series with a dash-pot and the Kelvin-Voigt (KV) model, where a spring and a 
dash-pot are connected in parallel mm- These two models provide the simplest 
representation of a visco-elastic material. The Maxwell model behaves like a liquid 
on long time scales, while the KV model represents a material behaving like a solid 
over long times. 

An inhnite number of elements connected in hierarchical patterns lead to power- 
law relations - the hallmark of fractal systems as discussed earlier by Heymans(2003). 
Fractional or generalized calculus is an alternative elegant method of tackling such 

systems mm- 

We try to construct an appropriate model keeping in mind the salient points of 
the results of our experiments. The results show the following interesting features (a) 
there is a hnite difference between the long-time strain for the single-step loading and 
two-step loading experiments. The difference increases with the time lag r between 
the two loading steps, (b) immediately after loading, the initial strain shows slight 
oscillations as reported earlier [6],(c) the strain tends asymptotically towards a more 
or less constant value in all three experiments, showing a solid-like nature of the 


material at long times. 


Several visco-elastic models are tried out to see whether these observations can 


be reproduced qualitatively. A 3-element solid-like model (the Boltzman model [15]) 
is found to give the most satisfactory £t to the experimental data provided certain 
features are incorporated into it. These are the following - a non-Newtonian nature 
is attributed to the dash-pot fluid through a non-integral time derivative in the 
viscosity equation; a change in the effective elastic properties of the material, similar 
to work hardening is introduced into the model. We allow the change in elastic 
modulus to increase with the time lag before the next loading step. 

we describe below three different models and compare their results with our 
experiments. The models are - (i) the generalized Kelvin-Voigt model, (ii) a three 
parameter fluid model and (hi) a three parameter solid model (Boltzman model). 
Fig [3] shows schematic diagrams of the models (i), (ii) and two equivalent versions 
of model (hi). 

4.1 (i) The generalized Kelvin-Voigt model 

It has been shown previously [6] that a Kelvin-Voigt element with a non-Newtonian 
viscous contribution as shown in £gure([3li) can reproduce the time development of 
strain shown by Experiment 1, with an initially oscillating strain which saturates 
asymptotically. Hence we begin our analysis of the present data using this model. 

In this case, the total stress a given below, is the sum of two terms, the hrst for 
a non-Newtonian dash-pot and the second for the spring. 


a = l3tc^- -h Ee 


(3) 
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Figure 3: Schematic diagrams of (a) the generalized Kelvin-Voigt model, (b) gen¬ 
eralized three parameter fluid model, (c) generalized three parameter solid (Zener) 
model, (d) generalized three parameter solid model(Boltzman model) 











































Here e is the strain and E is the elastic modulus. A parameter f3 with dimension of 


stress and a characteristic time characterizes the effective viscosity of the non- 
Newtonian fluid. For the Newtonian case q would be 1 and I3tc = rj, the coefficient 
of viscosity. 

Solving according to the procedure outlined in Dutta Choudhury et al.(2012) we 
arrive at the following expression for the strain. 




where, IE(— 2 ) is the one parameter Mittag-Leffler function dehned by 


Ei(-^) = e-" 



(5) 


Introduction of a non-integer order derivative q changes the dimensions of the 
physical quantities in the strain equation, so it is convenient to use non-dimensional 
parameters. The time has been scaled by a characteristic time tf. and the strain 
represented as a function of the scaled time e{t/tc)- The other ratios occurring in 
equation(j4j) - ^ and j are dimensionless. 

For step loading a stress (Tq is applied by placing the weight Wi at t = 0'*' which 
produces strain e = a^Jit). If at t = E some additional load IF 2 is applied there is 
a sudden change in the strain e(t)[IT]. The total strain may be written as 


e{t) = croJ{t) + AaJ{t — ti) 


( 6 ) 


where J{t) is the creep compliance at time t. 





Hence the final strain equation for two-step loading is, 


e{t/tc) = 


a 




a 

E 




(7) 


It is to be noted here that the load on the sample due to the weight W acts 
vertically, but the sample spreads out laterally under a shearing stress, so a in the 
equations H] and [7] is not numerically equal to the corresponding weight W, but we 
may assume it to be proportional to W. The proportionality constant is absorbed in 
the denominator E which is also proportional but not equal to the elastic modulus of 
the material. This argument is discussed in detail in Dutta Choudhury et ah (2012). 
Fig dUi) shows results using this model with q = 1.2. The initial oscillations observed 
when the load is placed on a drop of gelatinized potato starch is a signature of the 
non-Newtonian viscosity of the fluid |6]. 

Theoretical results for the strain calculated from the non-linear K-V model are 
depicted in Fig.(jl^). A problem with numerical methods in fractional calculus 
formalism using the Mittag-Leffler function is that beyond 30 time-steps, rounding 
off errors accumulate and large fluctuations are observed, as discussed by Dutta 
Choudhury et al.(2012) and Heymans and Bauwens (1994). Up to 20 time-steps 
results are reliable and realistic. The time-step is scaled accordingly, so that the 
dimensionless time remains within this range. 

A significant feature of our observed squeeze flow results is the ‘memory’ retained 
by the material of its loading history. We can see that the KV model does not re¬ 
produce this aspect. It cannot reproduce the dynamic visco-elasticity measurements 







Figure 4: The three hgures from left to right, show the calculated strains for the KV 
model, the 3-parameter fluid model and the 3-parameter solid (Boltzman) model 
respectively. The black solid line represents one-step loading and dotted (red) and 
dashed (blue) two-step loading with r = 5.3s and r = 12.3s respectively. 

either as evident from fig. ([5^), so obviously a more complex model is needed. We 
try a 3-parameter liquid model, which can lead naturally to a difference between 
one time loading and step loading in the long-time limit. 

4.2 (ii) Three parameter fluid model 

When a dash-pot is attached to the Kelvin-Voigt model in series, the resulting 
system represents a three-element fluid model and hence allows creep [lU. A linear 
version of the dashpots will not lead to the observed oscillations in the strain. To 
introduce non-linearity in the viscous part of 3 parameter fluid model, we replace 
the linear dash-pots by fractional order dash-pots, the outer dash-pot incorporates 
a fluid with the fractional order qi for the time derivative and the dash-pot in the 
K-V model is substituted by a q 2 order dash-pot. So the fractional order constitutive 





















Figure 5: The three hgures from left to right, show the calculated results for the 
storage and loss moduli for the KV model, the 3-parameter fluid model and the 
3-parameter solid (Boltzman) model respectively. The black solid line represents G' 
and the red dash-dot line represents G”. 

equation reduces to. 


WQl ( 1<12 

1 + {NE)— + {ME) 


dti^ 




a{t/tc) = 


fill rlQi+q2 


( 8 ) 


where,/9i and /32 are the apparent viscosities of outer dash-pot and the dash-pot 
in K-V model of hgure ([3 ]d) and E is the elastic modulus of the spring part of the 
model. 

Taking it’s Laplace transform and solving we get the strain-time relation as 

e(t/tc) = ao//3i + ao/l32E2{q2, a) (9) 

7(gi + 1) 


where E 2 (g 2 ,<a) is the second order Mittag-Leffier function dehned as , 
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The equation for the hnal strain after total loading, 


e{t/tc) = (Jo/13 


''1 


S(?i +1) 


+ ao//32E2{q2, a) + ao//3 


S(?i +1) 


+ ao//32El2{q2,(y) (H) 


Here 


El2{q2,cy) = ^ 

k=0 


r[(/c +1)^2 + 1 ] 


( 12 ) 


i.e, t in (ITOi) is replaced by {t — ti). 

For step loading we again follow the treatment described earlier to get the plot of 
e{t/tc) against t as shown in £gure(l4]3). The graph shows that the system remembers 
the loading history of the material i.e. in the asymptotic long-time limit there 
remains a gap between the strain curves for one-time loading and step-loading. 
However, the creep introduced by this model is not observed in the experiments, 
rather the strain saturates to a constant value for both one-time loading and step¬ 
loading. As these results do not ht our experiments satisfactorily, we revert to a 
solid-like model with some modihcations as described below. 


4.3 (iii) Three parameter solid (Boltzman) model with non¬ 
linear viscous elements 

The third model we consider is the Boltzman model [15] a 3-parameter visco-elastic 


model where the Kelvin-Voigt solid and a spring are connected in series as shown in 






Model 

For Full 

1st step-loading 

2nd step-loading 

2nd step-loading 

parameter 

loading 

r=0s 

at r=5.3s (A) 

at r=12.3s (B) 

a 

1.56 

0.78 

0.78 

0.78 

Q 

1.2 

1.2 

1.2 

1.2 

ai 

0.6 

0.6 

0.6 

0.6 

hi 

8 

8 

8 

8 

m 

1.13 

1.13 

1.30 

2.78 

a 

1.6 

1.6 

1.6 

1.6 


Table 1; The parameters for the 3-parameter Boltzman model under different load¬ 
ing sequences. 

fig- (EJl). This model has been shown to be equivalent to the Zener model [fig.([3t:)]. 
Such models have been used to describe visco-elastic behaviour of biological mate¬ 
rials such as Cartilage and white blood cell membrane [T5] . 

The Boltzman model shows visco-elasticity with long-time solid like behaviour. 
This model can be shown to be equivalent to the Zener model, which also consists 
of a dashpot and two springs ( see fig.([3]^c and d))). Our trials with the previous 
models lead us to the conclusion that we cannot model the memory of loading 
history using a solid-like model, assuming a constant set of model parameters. The 
phenomenon of strain-hardening or work-hardening, where visco-elastic properties 
change on loading is well known We now introduce this aspect into our 

visco-elastic model. The primary idea is that during the second step for two-step 
loading, the material has already been under stress for time r and this causes a 













change in its elastic modulus, the change being proportional to the time interval. 
So, while the one-time loading and first step of the two-step loading are under 
identical conditions, a change in material properties has to be considered during 
the second loading step. The viscous element - the dash-pot is replaced by its 
intermediate Scott-Blair fractional system of order q |1]. This concept is similar to 
the sectioi im The constitutive equation for the fractional Boltzman model is. 


[1 + = ['^ + 


(13) 


with ai = m = and h = 




Eq+Ei ’ Eq+Ei Eq+E\ 

the various elements are shown in the figures ([St;) and ([3li). A parameter jS 
characterizes the effective viscosity of the non-Newtonian fluid and when g is 1, 
fdtc = r], the coefficient of viscosity. Proceeding as in section (I4.ip . we get the 
solution of (1T3|) . 


e{t tc) = ao— + (To 

0 i \m 




The material functions are. 


Jit/tc) — jg + Jl 


1 - Eg - (to)' 


(15) 


•^9 “ If’ ~ If] ’ ^ ~ ^ contains the relaxation time tc which is 

taken as 1 here. For getting the total strain after step-loading we follow equation 


and finally reach the solution as. 


e(l/lc) ^ cToJ(l/lc) + ActJ[(1 - 


( 16 ) 










We assume that work hardening or lethargy develops in the material in the second 
step during step-loading. This is taken into account by assigning a stiffer elasticity 
to the sample through a modihed value of m, which contains the combination of Eq 
and El. The best-£t parameters for model (hi) are given in table fl4.3p and calculated 
curves are compared with experimental results in hg.([T]). 


4.4 Calculation of dynamic moduli G'{u) and G''{u) 


Expressions for G* = G'{uj) -|- iG''{uj) for the three models described above are as 
follows: 


4.4.1 Fractional Kelvin-Voigt model 


We have calculated complex relaxation modulus from the Kelvin-Voigt fractional- 
model. The form of this modulus is: 


G* = {E + cos -f- i{Puj'^sin^) 


(17) 


Here G' = E + cos ^ and G” = (duj^sin^ are the elastic modulus and loss 
modulus respectively. 

4.4.2 3-parameter non-linear fluid model 

The second model i.e. the 3-parameter non-linear fluid model reveals complex fea¬ 
tures of elastic modulus and loss modulus as 


, B cos ^(o^ -|- -|- 200;"' cos ^) -|- Gu‘^^{a + cos ^) 

® =- D - 


(18) 


and 



(19) 


,, sin ^ + CK^ sin ^ + 2au'^ sin ^ cos + Cci;"'+^ sin ^ 

^ “ D 

where n = p+q, B = C = and D = Pfoj^P+P^uji{oji+a‘^+2a cos f) + 
2/9 i/32(cj”' cos + acj^ cos 


4.4.3 3 parameter Boltzman solid model 


Complex relaxation modulus after taking the Laplace transform of constitutive equa- 
tion firs]) for 3 parameter Boltzman model, 

= ‘t>{s) 1™= fg 

gives 


G*(a;) = G'{u)+iG"{u) 


where 


G'iu) 


m + (aim + hi)uj^ cos qn/2 + 
1 + 2aiu‘i cos qn/2 + 


( 20 ) 


and 


G''{u) = 


ipi — mai)uj^ sm. qn/2 


( 21 ) 


1 + 2aiuj'i cosg7r/2 + 

G' and G" are plotted in Fig.(j2]) for the three models. G" increases linearly in 
the hrst two models (a) and (b), and shows a peak in the 3-parameter solid model. 
G' is a non-increasing function of frequency in the hrst two models, but exhibits a 


complex variation in the third model with a generally increasing trend. 






5 Discussion: comparing theory and experiment 


Let us compare the performance of the three different models presented in repro¬ 
ducing experimental results. 

The extended KV model |3^], with non-linear viscosity, successfully explained 
the oscillating strain for single step loading in the earlier work |6]. However, being a 
solid like model it cannot lead naturally to a long-time difference between single-step 
and multi-step loading, under identical total load m- The dynamic viscosity trends 
are also totally different from the experimentally obtained curves. So obviously a 
more complex model is needed to reproduce the new results. 

The three-parameter fluid model ([SJd) can reproduce naturally the asymptotic 
difference between single-step and multi-step loading. But, being essentially a liquid¬ 
like model, it shows creep. The strain continues to increase approximately linearly at 
long times, this is not in conformity with experiments. The experimental storage and 
loss moduli also show a dominance of the elastic i.e. solid-like behaviour compared 
to the dissipative, i.e. loss modulus. So, even with three elements and non-linear 
viscosity the three-parameter fluid model fails to reproduce the features seen in the 
real system. 

It is now obvious that to simultaneously exhibit a memory of loading history 
and long-term saturation of strain, a solid-like model with additional modification 
is necessary. A lethargy effect or strain hardening has been observed in polymeric 
as well as crystalline materials (TB] and is quite likely to occur in a complex material 
like a starch gel. The non-linear viscosity in the model leads to the oscillatory 


variation seen after initial application of the load. Introducing these in the three- 
parameter Boltzman solid model (d) we get strain variation curves (represented by 
lines) which can reproduce the experimental results (denoted by the symbols) quite 
well as shown in fig ([T]). The dynamic response of the system under oscillatory strain 
has been calculated for the three models and compared with experiments. Figures 
dSD show that experimental variations of storage and loss moduli do not match any 
of the calculated results exactly. However, models (a) and (b) yield very simple 
curves for G' and G” with cu where there is no indication of any increasing trend 
in G^(omega). Model (d) i.e. the modified three-parameter solid model shows a 
general rising trend for G'(omega), particularly noticeable in a certain frequency 
range. Though this model does not lead to the more-or-less linear rise seen in the 
experiment, it is closer than the other models. G” does not exhibit a peak in the 
experiment, as seen in model (d) results, but it is possible that the frequency range 
accessible in the experiment is inadequate for this. Evidently, more modifications 
of the model are needed, for faithful quantitative representation of dynamic visco¬ 
elasticity behaviour. For the present we may consider model (d) a good starting 
point. Other reports on potato starch and other starches show a behaviour similar 
to our experimental results [TO] . 

6 Conclusions 

The interesting visco-elastic behaviour of starch gels, particularly potato starch is 
manifested in several very simply performed experiments. Micrographs show potato 


starch to consist of rounded cellular structures dZllIB] which are known to contain 
amylose and amylo-pectin. On heating, the long polymer chains of amylose leach 
out of the cells,causing a thickening termed gelatinization. The highly viscous fluid 
produced exhibits non-linear viscoelasticity and other forms of instability. For ex¬ 
ample, finger-like undulations develop at the periphery of a squeezed drop! [T9] 1 and 
there are oscillations in the area of contact of the drop and substrate. Potato starch 
gel also acts as a host for crystal growth and a drying drop of the gel with added 
NaCl can produce salt crystals with a wide range of morphologies [20] . 

These results are important in food industry ra. related biological systems pi] 
as well as from the general point of view of visco-elasticity. Considering the host 
of interesting phenomena observed in such very simple experiments, studies of the 
type reported here promise to be very useful. 
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